Statistics of S-matrix poles for chaotic systems with broken time reversal invariance: 

a conjecture 
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In the framework of a random matrix description of chaotic quantum scattering the positions of 
S— matrix poles are given by complex eigenvalues Zi of an effective non-Hermitian random-matrix 
Hamiltonian. We put forward a conjecture on statistics of Zi for systems with broken time-reversal 
invariance and verify that it allows to reproduce statistical characteristics of Wigner time delays 
known from independent calculations. We analyze the ensuing two-point statistical measures as e.g. 
spectral form factor and the number variance. In addition we find the density of complex eigenvalues 
of real asymmetric matrices generalizing the recent result by Efetov [ jTj . 



PACS numbers: 05.45.+b 

One of the basic concepts in chaotic quantum scatter- 
ing is the notion of resonances, representing long-lived 
intermediate states to which bound states of a "closed" 
system are converted due to coupling to continua. On a 
formal level resonances show up as poles of the M x M 
scattering matrix S a b(E) occuring at complex energies 
£k = Ek — (i/2)Tk, where Ek is called the position and 
Tk the widths of the corresponding resonance, and M is 
the number of channels open in a given interval of en- 
ergies. Recently, advances in computational techniques 
made available resonance patterns of high accuracy to 
be obtained for realistic models of atomic and molecular 
chaotic systems 0. 

As is well-known, universal statistical properties of 
bound states in the regime of quantum chaos have their 
pattern in statistics of real eigenvalues of large random 
matrices The methods to adjust random matrix de- 
scription to the case of resonance scattering in open quan- 
tum systems are very well-known since the pioneering 
work by the Heidelberg group . The method proved to 
be very fruitful and allowed to calculate different univer- 
sal characteristics of chaotic scattering, see papers [|]|| 
for a thorough discussion of recent developments. 

In the framework of this approach the S-matrix poles 
(resonances) are just the complex eigenvalues of an ef- 
fective random matrix Hamiltonian Ti. e ff — H—iT. Here 
H is a large sclf-adjoint N x N matrix of appropriate 
symmetry serving to describe the statistical properties 
of the closed counterpart of the scattering system under 
consideration. The NxN matrix T is to describe a possi- 
bility of transitions from the states described by H to the 
outer world via M open channels and is related to the 
NxM matrix W of transition amplitudes as r = 7rWW''. 
Such a form is actually dictated by the requirement of 
S— matrix unitarity and ensures that all the S— matrix 
poles are in the lower half-plane of complex energies as 
required by causality. 

In spite of quite substantial analytical and nu- 

merical work on properties of Ti. e f / our actual knowl- 



edge of resonance statistics is rather limited. Apart from 
the simplest perturbative Porter-Thomas treatment p]] 
( as well as its more advanced variants [O) the results 
available on statistics of S— matrix poles are: (i) the 
joint probability density of all resonances for the system 
with one open channel and unbroken time-reversal invari- 
ance (TRI) 0, (ii) the mean density of S-matrix poles 
in the complex plane for large number of open channels 
M ~ N ^> 1 ||, (iii) the mean density of S-matrix poles 
for arbitrary M <C -/V for the case of broken TRI f§ . In 
particular, no information about two-point correlations 
between different poles is available to our best knowl- 
edge. 

The situation improves drastically if one replaces the 
physically motivated matrix T introduced above by an 
(unphysical) NxN Hermitian matrix A with indepen- 
dent, Gaussian distributed elements and considers H to 
be taken from the Gaussian Unitary Ensemble (GUE), 
thus restricting attention to the systems with broken 
TRI. When the variance of both H and A coincide, 
such an ensemble was studied long ago by Ginibre [ p"2| , 
who managed to find all the correlation functions of the 
eigenvalues in the complex plane. They turned out to 
be quite different from those known for the self-adjoint 
Gaussian random matrices with real eigenvalues, studied 
by Wigner, Dyson, Mehta and others |Q. 

At the same time it is clear that reducing the vari- 
ance of A as compared to that of H drives the ensem- 
ble towards GUE. The existence of a nontrivial regime 
of weak non-Hermiticity was recognized in our preced- 
ing works This 

regime occurs when TrA ~ 
jjTrH 2 in the limit of large N. Under this condition 
the imaginary part Yk of a typical complex eigenvalue 
Zk = %k — iYk is comparable with the mean separation 
A = (xk — Xfc+j) ~ jf between neighboring eigenvalues 
along the real axis. Exploiting the method of orthogo- 
nal polynomials we demonstrated jlij that all statistical 
properties of H. = H — iA in this regime can be described 
in terms of a kernel K(Zi,Z2) depending on two com- 



1 



plex coordinates Z\ t % = x±ui/2N— iy-i^/N. In particular, 
the mean eig envalue density p(Z) = ( J2k ^\ z ~ Z h)) 
is given by p{Z) = K(Z, Z) and the two-point cluster 
function y 2 {Z\, Z 2 ) defined via the relation: 

(p(z 1 )p(z 2 )) c = (p(z))8^\z 1 -z 2 )-y 2 (z 1 ,z 2 ) (i) 

is given by y 2 (Z u Z 2 ) = \K{Z U Z 2 )\ 2 . 

Unfortunately, such a detailed information is of no 
direct use for the physically motivated case of chaotic 
scattering. Nevertheless, the insights provided by the 
Gaussian case combined with the existent knowledge of 
the resonance statistics allowed us to put forward a well- 
grounded conjecture about statistics of complex eigenval- 
ues of weakly non-Hcrmitian ensemble of the type H—iY 
for any given Hermitian matrix T. This issue is the con- 
tent of the main part of the present paper. 

Before formulating the conjecture, let us recall that the 
GUE ensemble is an invariant one, i.e. its statistics is the 
same independently on the basis chosen. Therefore, one 
always can go to the eigenbasis of the matrix F and con- 
sider it to be diagonal with eigenvalues ji, i = 1,...,N. 
In what follows we find it convenient to characterize the 
matrix F by the following function: 



/r(z)=5>(l 
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^ij(7i and v{x) 



where gi 

stands for the Wigner semicircular density of real eigen- 
values of the Hermitian part H of the matrices Ti.. 

Now we put forward the following conjecture: 
Suppose the function fr(z) defined above has a finite 
limit when iV— >oo. Then the statistics of eigenvalues Zi 
of the corresponding almost Hermitian ensemble H e f / in 
the limit N 3> 1 is completely determined by the kernel 



N 2 



■kv sc {x) 



K{Z 1 ,Z 2 )^—e 1 ^- y ^ I due- u{vi+V2)+iuuJ+fT{u) (3) 



-7Tf sc (x) 
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In particular, the mean eigenvalue density p(Z) = 
( J2k Si2) ( z ~ z k)} is given by p(Z)=K(Z,Z) and the 
two-point cluster function is y 2 (Zi, Z 2 ) = \K(Z\, Z 2 )\ 2 . 

Let us now systematically veryfy the compatibility 
of our conjecture with the known properties of almost- 
Hermitian matrices of various types. 

The simplest test is to make sure that for a Gaussian T 
such that TrT 2 ~ -^TrH 2 we are back to results proved 
in JT3] . Indeed, for this case a typical eigenvalue 7, ~ 
N-^T 2 <^ i_ 5 hence g~ x « 2iTv{x) r )i <C 1 and in the limit 
of large N one can expand f(z) in a series. The first term 
(proportional to z) vanishes in the limit N^oo because 



of the symmetry of the distribution of eigenvalues of f 
around zero. Thus, the leading term turns out to be 
proportional to z 2 . The corresponding Gaussian integrals 
over k\^ 2 in Eq.(|j|) can be performed exactly, the resulting 
kernel reproducing that found in |l4[ . 

Let us now consider the mean eigenvalue density 
p(Z) = K(Z,Z). The "physical" case L = nWW^ with a 
finite number M of open channels (i.e, with a finite num- 
ber M of positive non-zero eigenvalues Ji, i = 1, --,M) 
was considered earlier in Q. The result coincides with 
that following from Eq.(||). 

Actually, one can easily adopt the methods used in || 
to satisfy oneself that the validity of the corresponding 
expression is not restricted by the case of positive 7$, but 
rather extends to an arbitrary set of eigenvalues. This 
fact provides a proof of our conjecture on the level of the 
mean eigenvalue density. 

Let us now show that our conjecture survives a much 
more stringent test on the level of two-point correlations. 
For this purpose let us invoke the notion of the so-called 
(energy dependent) Wigner time delay defined in terms 
of resonance positions Ef. and widths as (see || for 
more details): 
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(4) 



Using this expression it is easy to relate the correla- 
tion function (tw {E\)tw {E 2 )) c of the Wigner time de- 
lays at two different energies Et t 2=E±Sl/2 to the two- 
point correlation function (p(Zi)p(Z 2 )) c of the densities 
of S— matrix poles in the complex plane Z — x— iY . Con- 
sidering the energy difference E\ — E 2 = Q, to be com- 
parable with the mean level spacing A = l/v(E)N and 
exploiting the fact that both the mean density (p(x,Y)) 
and the cluster function y 2 (xi, x 2 , Y\, Y 2 ) change with 
x— (xi+x 2 ) /2 on a scale much larger than A, one can 
perform the a;— integration explicitly. After this it is 
convenient to pass to the scaled variables: 



T W 



nY 



TTUJ 
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with ijj = x\— x 2l to rescale the cluster function and the 
density as follows: 

PE{y) = — {p{E 1 Y)) ] y 2 {E,uj Xl y ll y 2 )^^y 2 {E 1 u Jl Y ll Y 2 ) 

and to make a Fourier transformation with respect to 
Co. These manipulations allowed us to write down the 
relation we were looking for in quite a compact form: 

00 00 
\JdCoe lM (r w {E+^Co)f w {E-^Co}j = Jdye~ 2ty p E (y) 



du x Jd yi Jdy 2 e- t{ y^+^y 2 {E,uj x , yi ,y 2 ) (5) 
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Such a relation between the correlation function of 
time delays and two-point cluster function provides a 
possibility of the most non-trivial test of our conjecture. 
Indeed, both the left-hand side and the first term in the 
right-hand side are known from independent calculations 
||, which allows us to rewrite the Eq.(||) for t > as: 



oo poo poo 
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z J— min(t,l) i=1 



9i + A 



g t + X + t 



(6) 



Although this relation does not provide a possibility to 
extract the cluster function in a unique way, it is easy to 
satisfy oneself that a direct subtitution of y 2 (Zi, Z 2 ) = 
\K(Zi, Z 2 )\ 2 , with the kernel taken from Eq.(||) into the 
left-hand side of Eq.(^|) produces after rescaling and in- 
tegration exactly the right-hand side of Eq. . This fact 
provides the strongest support to our conjecture. 

Having at our disposal the conjectured form of the clus- 
ter function y(Z\, Z 2 ) it is interesting to calculate other 
related quantities frequently used in applications, such as 
the spectral formfactor and the number variance. For the 
Gaussian case such a calculation was performed in [fTi) . 
For the physically interesting case of resonance scatter- 
ing the formfactor can be obtained along the same lines 
and is equal to: 



b(q u q 2 ,k) = I dw x jdy 1 jdy 2 e 2 ^^ + ^ +k ^y 2 {Z 1 ,Z 2 ) 

(g t is + v) 2 — k 2 



-oo 
{v-\k\) 
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1 {giV+v-iq 1 )(g i v+v-iq 2 ) 



(7) 



with v = v{X). The variance of the number of resonances 
in a strip < RcZ < L; — oo < ImZ < oo of the width 
L = L X A comparable with A can be expressed in terms of 
&(0, 0, fc), see Q. In particular, for the case of equivalent 
scattering channels: gi—g, i = 1, M it is given by: 

If 1 

^2(L X ) = L x ^ / dkk 2 sin 2 (irkL x ) 

7T JO 



\ / dv[l-k 2 {g + vy 2 ] 

-(l-fc) 



M 



(8) 



where we have used ttv(0) = 1. Let us discuss for sim- 
plicity its typical features for the symplest case of only 
one open channel M = 1. We are actually interested in 
deviations of the number variance from its value known 
for Hermitian GUE matrices. One finds: 



6Y,{L x ) = V 2 {L x )-Y,f UE \L x ) 
2 f 1 

= — / dk sin 2 (7rfcl/ x ) 



(9) 



(l-fc) 



9 2 ~ (1 - kf 



Usually, one is interested in the behaviour of the num- 
ber variance for L x 3> 1. Then, for any g > 1 one 
finds that the difference 8Yi(L x ) tends to a constant value 
In [g 2 /(g 2 — !)]■ The so-called "perfect coupling" case 
g = l is known for being specific: e.g. the resonance den- 
sity shows at g = 1 the powerlaw behaviour || . On the 
level of number variance such a specificity is reflected in 
a logarithmically growing difference: 5Y,(L x )cx\iiL x . 

Finally, the knowledge of the cluster function allows 
one to find the small-distance behaviour of the so-called 
nearest- neighbour distance distribution p(Zq, S<A), see 
e.g. Q. As expected 0, the leading term for S — > 
turns out to be always cubic: p(Z 0} S) ~ S 3 , as long as 
the system is open: <?~ 1. However, it is interesting to 
note that for very weakly open systems: (?3>1 the cubic 
law is modified in such a way, that in a parametrically 
large region g^ 1 -C 5/ A < ImZo/A the behavour becomes 
p(Z Q , S) - S 5 / 2 (cf. JH). Unfortunately, the density of 
complex eigenvalues is exponentially small under the con- 
dition ImZo/A3> g^ 1 , the fact hindering a possibility to 
observe such an unusual regime numerically. 

Let us now turn our attention to another class of ran- 
dom matrices with complex eigenvalues which attracted 
much attention recently. Namely, we consider the en- 
semble of weakly asymmetric matrices with real elements. 
Such matrices can always be presented in the form H+A, 
with H being real symmetric (hence, taken from GOE) 
and A being a real antisymmetric: Ay = — Aji such 
that NTtA 2 - TtH 2 in the limit of large N. The case 
of matrices A with independent, identically distributed 
Gaussian entries was studied by various authors and by 
different methods jl7|,[l6 15 1. In particular, the follow- 
ing unusual property was detected numerically in Jl6| ] 
and proved analytically in p5|Jl7| : the finite fraction of 
eigenvalues stays on the real axis even for A^O. This 
fact should be contrasted with the corresponding prop- 
erty of earlier discussed weakly non-Hermitian matrices 
whose eigenvalues with probability unity have a finite 
imaginary part as long as f ^ 0. More recently, the 
interest to the ensemble of slightly asymmetric real ma- 
trices arose after the work by Efetov |p^| , who discov- 
ered its relation to an interesting problem of motion of 
vortices in disordered superconductors with columnar de- 
fects fl8[| . Efetov calculated explicitly the mean density 
of eigenvalues for the Gaussian A. Shortly after Halasz et 
al discovered that Efetov's result describes also the 
density of real eigenvalues of some matrices appearing 
in a random matrix approach to the problem of sponta- 
neous breaking of chiral symmetry in QCD. An interest- 
ing feature was that the perturbation considered in M] 
forcing the eigenvalues to leave the real line was not at 
all random. Translating these results to the ensemble 
of random real non-symmetric matrices it is natural to 
expect that for antisymmetric perturbations of the type 
/il 
-fil 



A = 



with a constant /i being of the order 
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of 1/ y/N Efetov's formula should still provide the correct 
eigenvalue density. 

This fact motivated us to reconsider the problem of 
calculation of the mean eigenvalue density for a general 
fixed real antisymmetric matrix A as it is done above for 
the case of almost-Hermitian matrices. Invoking again 
the arguments of the rotational invariance, it is enough 
to consider the matrices A of the following structure: 
A = diag(Ai, An), with each block Ai being 2x2 ma- 
fj,i 
-Mi 

symmetric A can be reduced to such a form by an orthog- 
onal rotation. The density of complex eigenvalues for the 
matrix H+A can be found by a straightforward modifica- 
tion of the calculation presented in [JL7| . Introducing the 
scaled variable y = iw(X)Nlva.Z and rescaling the den- 



trix of the form A = 



, since an arbitrary anti- 



sity correspondingly: px(y) 
finds: 



(p(Z))/(NMX))\ 



p x ( y )=5{y) due ilM^u)+f A (^u)] ( 1Q ) 
Jo 
1 f 1 

+ _ / duusmh(\y\u)e? [M ™ u)+fA{ - 7nyu)] 
2tt Jo 

oo />oo 

ds / rffce~ ifes ~5[/- 4 ( iTriyfc )+/^(- i7r,yfe )] 

where the function /a(z) is given again by Eq.(^) with 
7i replaced by /Xj. 

From the derived expression one immediately infers 
that if a typical ^ is of the order of iV -1 / 2 , the func- 
tion can be expanded up to a first nonvanishing 
order such that Ja{u) + fj^(—u) ex TiA 2 u 2 and the cor- 
responding expression coincides with that found by Efe- 
tov. As such, Efetov's formula is indeed applicable also 
for constant matrices A of the type described above. 

We see that the most striking qualitative features of the 
Efetov's formula: i)a nonvanishing density of real eigen- 
values and ii) a linear increase with \y\ of the probability 
density to have a finite imaginary part -persist for any 
antisymmetric perturbation A. The strongest quantita- 
tive deviation from Efetov's result occurs in the case of 
finite-rank perturbations A such that pn — for i > M. 
In particular, one can show that if at least one of the 
quantities gi = ■^^{l 1 i J r l 1 ^ 1 )- 1S equal to unity, the mean 
density decays asymptotically as p(y 3> 1) oc y~~ 2 . Such 
a slow power law decay should be contrasted with the 
Gaussian case when one always has a very sharp cutoff 
of the density for large enough y. For the general case of a 
finite-rank antisymmetric perturbation such that gi =/= 1 
the density is cut exponentially at y ~ (gi — 1)~ . 

In conclusion, we put forward a conjecture on statis- 
tics of S-matrix poles Z^ for systems with broken time- 
reversal invariance and verified that it is perfectly com- 
patible with the existent knowledge on quantum chaotic 
scattering. In particular, our conjecture allowed us to re- 
produce statistical characteristics of Wigner time delays 



known from independent calculations. We analyzed the 
ensuing two-point statistical measures as e.g. spectral 
form factor, number variance and small distance behav- 
ior of the nearest neighbor distance distribution p(Zq, S). 
In the final part of the paper we calculated the density of 
complex eigenvalues of an ensemble of real weakly asym- 
metric matrices. The expression obtained generalizes the 
recent result by Efetov jlTj to the case of an arbitrary 
antisymmetric perturbation. 
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